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Abstract 

We probe the dynamic strength of multiple biotin-streptavidin adhesion bonds under linear load- 
ing using the biomembrane force probe setup for dynamic force spectroscopy. Measured rupture 
force histograms are compared to results from a master equation model for the stochastic dynamics 
of bond rupture under load. This allows us to extract the distribution of the number of initially 
closed bonds. We also extract the molecular parameters of the adhesion bonds, in good agreement 
with earlier results from single bond experiments. Our analysis shows that the peaks in the mea- 
sured histograms are not simple multiples of the single bond values, but follow from a superposition 
procedure which generates different peak positions. 
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I. INTRODUCTION 



Cell adhesion in multicellular organisms is mediated by highly specific, but weak receptor- 
li Grand bonds QB.Q . Typical bond energies are in the range of several /c#T so that bond 
lifetimes are finite due to thermal activation from the environment. Since cellular adhesion 
sites usually operate under mechanical load, different experimental techniques have been 
developed to study rupture of single biomolecular bonds under an applied load. Motivated 
by atomic force microscope (AFM) measurements on biotin-streptavidin [4], it has been 
shown that the force at which a bond breaks under an applied load is a stochastic variable 
with a probability distribution which depends on the loading protocol and which can be 
explained by Kramers theory for the escape of an overdamped particle from a potential well 
over a sharp transition state barrier (3, tsj] . This insight led to the new field of dynamic 
force spectroscopy (DFS) For the biotin-streptavidin system, the concept of dynamic 
bond strength has been confirmed by experiments with the biomembrane force probe (BFP) 
8j. If force F increases linearly in time t with loading rate m, i.e. if F — mt, a rupture 
force distribution results in which the most frequent rupture force F* is proportional to the 
logarithm of loading rate (the same result has been found for the average rupture force j^]). 
Dynamic force spectra use this relation to chart the binding landscape: a straight line in 
the plot of F* versus lnm is characteristic for a single dominant barrier limiting escape. 
The slope is determined by the distance of the barrier from the ground state (the reactive 
compliance Xb) and extrapolation to F* = yields the unstressed off-rate ko of the bond 
over this barrier. A series of energy barriers reveals itself as a sequence of straight lines with 
increasing slope, thus defining different regimes, each dominated by one transition state 
barrier. 

Single bond DFS experiments require a very low frequency of successful binding events to 
avoid multiple attachments. This renders single bond experiments time consuming and accu- 
mulation of rupture force histograms difficult. In the following, we describe DFS-experiments 
with the BFP on the biotin-streptavidin system at four different loading rates, in which we 
allow formation of multiple bonds. Experimentally, the streptavidin-biotin system is the 
best studied example for molecular bonding. In detail, it has been investigated with AFM 
4], BFP jg] and flow chambers [if]]. In the latter case, multiple bonds have been probed, but 
a persistent problem with flow chambers is that force is distributed in a heterogeneous way 
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over the different bonds, thus rendering a quantitative analysis difficult. In contrast, force 
distribution is expected to be significantly more homogeneous for the BFP as used here. 
The higher frequency of binding events introduced here makes acquisition of rupture force 
histograms considerably easier. However, in order to evaluate the histograms, a theoretical 
model is required which allows to extract the single bond properties from the multiple bond 
data. 

Deterministic models for the rupture of multiple, parallel bonds have been analysed before 
for constant [1] as well as linear loading [loj]. While deterministic models describe the 
average number of bonds and are most appropriate for large systems, stochastic models are 
required for small systems with a finite number of bonds. Here the simplest possible case 
is the irreversible rupture of a small number of equivalent bonds, which can be described in 



the mathematical framework of a one-step master equation |9j. Multiple bonds a. 
rebinding of broken bonds, which can be described in the same framework 11 



so allow 



121. It has 



13] . Similar 



also been applied before to the case of linear loading commonly used in DFS 
approaches have been also used to describe DFS on titin, where the different bonds may 
correspond to different hy drog en bonds within a single Ig27-domain [l^] or to the different 



Ig27-domains within titin 



In this paper, we use the framework of a one-step master equation to evaluate data from 
multiple bond DFS on the well-established streptavidin-biotin system. Because multiple 
bonds are now allowed in the experimental setup, the exact number of initial bonds in each 
experiment is a stochastic variable. By comparing numerical solutions of the master equation 
to experimental histograms, we can estimate the corresponding probability distribution. In 
this way, we show for the first time how the dissociation spectrum of multiple bond DFS 
results from the superposition of the contributions from different numbers of initially closed 
bonds. 



II. EXPERIMENTS 

To monitor rupture forces between multiple streptavidin-biotin bonds at different loading 
rates, we used the BFP-instrument as described previously [8] and depicted in fig. [H In brief, 
for use as force transducers, biotinylated red blood cells were pressurized into a spherical 
shape with a micropipette. The suction pressure sets the red cell membrane tension, and 
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FIG. 1: Videomicrograph of the biomembrane force probe (BFP). The red blood cell (on the right) 
acts as an elastic element, which transforms displacement of the pipette into force exerted on the 
binding site. The probe bead (attached to the red blood cell) and the test bead (on the left) are 
glass spheres (3 /xm in diameter) decorated with complementary receptor and ligand molecules, in 
this case streptavidin and biotin. 

thus the BFP spring constant, kp, in a tunable fashion at values as low as 0.1 — lpN/nm. 
Probe and test glass microspheres were decorated respectively with streptavidin and bi- 
otin following a common procedure, in which polyethylene glycol (PEG) crosslinkers (Nek- 



tar) were used to inhibit non-specific adhesion [171]. The streptavidin-coated probe bead 
was firmly attached to the apex of the biotinylated red blood cell capsule. The biotin- 
coated test bead was maintained in a second micropipette and manoeuvred to/from the 
BFP probe at controlled impingement forces (< 25 pN) and retraction speeds by precision 
piezo-translators. Bead positions were determined by on-line video processing at a sampling 
rate up to 180 frames/s and with a spatial resolution of about 6 nm. Validation of our 
set-up was achieved by measuring rupture forces of individual streptavidin-biotin bonds for 
loading rates between 5 and 50,000pN/s (data not shown, see 18j for details). The ob- 
tained dynamic force spectrum, which collects the most probable unbinding force for each 
loading rate, was in excellent agreement with previous reports js| . In particular, the derived 
energy landscape exhibits two transition barriers characterised by two reactive compliances, 
Xb,i = 0.14 nm and Xb.2 = 0.51 nm. 

As compared with single bond DFS, we mainly implemented two changes in the exper- 
imental procedure. First, while an attachment frequency of 1 per 10 touches is the usual 
criterion that ensures high probability of single bond formation, we increased the surface 
density of specific sites on the glass beads to reach an adhesion frequency of about 1 per 
3 touches. According to Poisson statistics, this increases the likelihood of picking multiple 
bonds. Second, and more important, since our goal was to quantitatively account for the 
distribution of all measured rupture forces, good accuracy for both low and high forces was 
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FIG. 2: Schematic representation of an adhesion cluster in the DFS experiment from fig.[U Between 
substrate and transducer Nt receptor-ligand bonds were established (here iVj = 5) of which i are 
still closed (here i = 3) while Nt — i have already ruptured (here Nt — i = 2). The load F applied 
to the transducer is assumed to be equally distributed over the i closed bonds, thus each of them 
feels the force Ff, = F/i. In our experiments, force is increased linearly in time with loading rate 
m, i.e. F = mt. 



required. Despite its advantages of tunability and softness, the red blood cell-based trans- 
ducer is only linear for extensions below ~ 0.5 /im 19|. In practice, this led us to set the 
BFP spring constant around lpN/nm and the range of explored loading rates between 100 
and 5000 pN/s. In consequence, we focused on the regime dominated by the outer barrier 
of the streptavidin-biotin energy landscape (with a reactive compliance Xb,2 — 0.51 nm). 

For biotin-streptavidin, rupture experiments with different experimental techniques in- 
dicated a dependence of unbinding forces on binding history, namely on the contact time 
that is available for binding before loading begins. In [20] this was assigned to the existence 
of an additional energy minimum (the absolute ground state) which is reached only after a 
large contact time and is more stable against rupture than the intermediate ground state 
that is reached first. In our experiments, we use the same contact protocol as in [8] with 
short contact times of 2 — 3 s, thus we expected the bonds to unbind from the same initial 
state probed before. 
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FIG. 3: Rupture force histograms for loading rates (a) m = 130, (b) 615, (c) 1250 and (d) 5000 pN/s 
which include 220, 170, 230 and 228 attachment events, respectively. The shaded bars show 
experimental results which were collected in bins of approximately 12 pN width. The histograms 
are normalised in respect to the total number of events and the bin width, so that they show the 
probability density of events. The solid lines are the simulation results fitted simultaneously to the 
complete data set (a-d). For the fitting procedure, the numerical results were discretized. These 
discretized curves are shown as dashed lines. The optimal fit parameters are ko = 0.02 Hz and 
F = 7.83 pN. 

III. MODEL 



Fig. [2] shows a schematic representation of the experimental setup from fig. [TJ Between 
two opposing surfaces N t receptor-ligand pairs are arranged in parallel. Each receptor can 
either be bound to its ligand (closed bond) or unbound (open bond). At a given time 
t there are i closed and N t — i open bonds. For the large loading rates that we use in 
our experiments, it is sufficient to consider irreversible bonds which cannot rebind after 
rupture [h], • Thus, the number Nt of bonds which is relevant for the rupture process 
is the number of closed bonds that has formed prior to loading. Rupture of a bond occurs 
stochastically through thermal activation. Following Bell [l| the off-rate for rupture of a 
bond increases exponentially with the force Fj, exerted to it, i.e. k ff = ko exp(Ff,/ F ). Here 
k is the unstressed off-rate at vanishing force and F = k B T/xb is the intrinsic force scale 
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for bond rupture, which is set by thermal energy fc#T = 4.1pNnm and reactive compliance 
x b . 

An essential part of the modelling is a reasonable assumption on how force is distributed 
over the closed bonds in a cluster. Here we assume that the force F = mt, which is de- 
termined by the displacement of the red blood cell, is independent of the number of closed 
bonds and is shared equally between them, i.e. F b = F/i. In our setup, this assumption is 
appropriate because the bonds are attached to a solid support of finite curvature. Because 
the binding times of tethered polymer bonds increase very rapidly with the separation of 
the opposing surfaces 21[ , our short contact times make formation of bonds at large separa- 
tions unlikely. Moreover biomolecular bonds carried by extended polymer tethers dissociate 
rapidly under force [22j. Therefore we expect the relevant bonds to form mainly within a 
narrow range of surface separations corresponding to the area of closest approach. Because 
the tethers that attach the ligands to the solid support have a contour length (~ 23.2 nm) 
which is small compared to the extension of the soft transducer, unbinding of a tether has 
a negligible effects on transducer extension and force. 

The number of closed bonds i in a cluster reduces with time from i = N t (initial, bound 
state) to i — (final, dissociated state). Because for shared loading all bonds are equivalent, 
the probabilities Pi(t) (i = . . . Nt) to find % closed bonds at time t completely characterise 
the system. Their time evolution follows the one-step master equation 

dPi m 
= r i+1 p i+1 - . (1) 

This equation states that i decreases through rupture of a closed bond with the rate = 
ik ff = iko exp(mt/iFo), where the linear factor i accounts for the fact that at any given 
time, each of the i remaining bonds can be the next to break. Once the master equation is 
solved, one can calculate the cluster dissociation rate D Nt (t) = po(t) = ripi(t) as a function 
of time, which can be converted to Dpf t (F) using the linear relation F = mt between time 
and force. 

Our model contains three unknown parameters, which can be determined through com- 
parison with experiment: the unstressed off-rate ko, the intrinsic force scale Fq (or, equiv- 
alently, the reactive compliance Xb) and the initial number of closed bonds Nt. In our 
experiments, the latter is itself a stochastic variable and cannot be controlled for every 
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FIG. 4: Superposition of the different Dpj t (F) (dashed lines) resulting in the final dissociation 
probability D{F) (solid line) fitting the experimental histogram for loading rate m = 1250 pN/s. 
The different DN t {F) are weighted with the distribution p a {Nt) shown as a histogram in the inset. 
The solid line with crosses is a least-mean-square fit to a Poisson distribution. 



rupture event. Instead we assume that N t follows a probability distribution 



Pa{N t ) — ocNtl ( a i H 1" "ivrO • 



(2) 



which is characterised by the set of iV t m coefficients a = {a\, . . . , ajy™}- is the maximal 
number of initial bonds considered and will be chosen in such a way that it does not affect 
the outcome of our parameter estimation. Because experimentally we only record data 
involving deformation of the transducer, events with N t = are not considered in the 
theoretical analysis. If closed bonds formed at a constant rate over a given time interval, 
p a (N t ) would be a Poisson distribution 



Px(N t 



A 7Vi e" A 



[l-e-*)N t \ 



(3) 



n 



as has been reported before for similar experimental setups 123. 1 . The Poisson distribution 
contains the single parameter A which determines the average number of broken bonds as 
(N t ) = \/(l-e-*). 

Experimental histograms are linear combinations of the distributions D^ t (F) for single 
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FIG. 5: Rupture force histogram with 230 data points and for loading rate m = 1250 pN/s 
simulated from the Poisson distribution from the inset of fig. HI The lines and the inset represent 
the same analysis as done before for the experimental data. 



N t in which the distribution eq. (jSJ) determines the coefficients: 



D(F)= ^p a (N t )D Nt (F). 

N t =l 



(4) 



For shared loading with a time- dependent force, analytical solutions for eq. ([I]) do not 
exist 13]. Therefore we use the Gillespie algorithm for exact stochastic simulations to 



generate rupture trajectories with the stochastic dynamics described by the master equation. 
For given values of k and F , distributions D Nt (F) of rupture force F are calculated for 
N t = 1 . . . N™ at the four loading rates. The distributions are superimposed via p a {N t ) as 
in eq. (j4j) and the relative weights are determined by a least-mean-square fit of D(F) to 
the complete set of experimental data. For this purpose, the simulation data are discretized 
in the same way as the experimental data and every bin is used as one data point. This 
procedure is iterated until values for fc an -d F are found that minimize the deviations 
between simulations and experimental histograms. 
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IV. RESULTS 



Fig. [3] shows experimentally obtained rupture force histograms including 220, 170, 230 
and 228 rupture events for the four different loading rates m = 130, 615, 1250 and 5000 pN/s, 
respectively. The histograms are compared to simulated rupture force distributions which 
are calculated using ~ 10 7 rupture trajectories at each value of N t and loading rate. Fitting 
simulations to measurements yielded the unstressed off-rate k = 0.02 Hz and the intrinsic 
force scale F = 7.83 pN, corresponding to the reactive compliance Xb = ksT/Fo = 0.52 nm. 
Thus our results extracted from multiple bond data compare very favorable with earlier 
results from single bond experiments, which gave Xb = 0.51 nm (a value for fc has not been 
given before). Closer inspection of fig. [3] shows that the simulated distributions approximate 
the positions and heights of the peaks in the histograms in a reasonable way. Note that the 
fits were not done separately to the different data sets for the different loading rates, but 
simultaneously to the complete set of experimental data. Thus the model is capable to fit 
well the whole range of loading rates probed. 

Fig. H] shows for the case m = 1250 pN/s how the different rupture force distributions 
D]y t (F) combine to give D(F). Each curve D^ t (F) displays a single maximum followed 
by a super-exponential decay at large forces. In general, their peaks cannot be identified 
with the peaks in the experimental histograms as it was done in previous analyses (e.g. in 
% [25|, [26J]). The best agreement is found for N t = 1, i.e. between the peak of Di(F) and 
the first experimental peak. In contrast, the positions of the peaks in D%(F) or D^(F) 
cannot be guessed from the histograms. In general, the peak positions of the D Nt (F) do not 
increase linearly with N t . Although this has been pointed out before [h]], a force quantum 
has occasionally been discussed in the literature. Our analysis shows that the situation is 
more complex and that accurate identification of the Dj^ t (F) from the simulation data is 
required to extract the single bond data from the multiple bond histograms. 

The extracted probability distribution p a (N t ) for the number of initial closed bonds N t 
is the histogram shown as inset in fig. HI Its average is (N t ) = 3.8. The largest Nt included 
was N™ = 10, because larger N™ did not make a difference for the estimate of the molecular 
parameters. The inset also shows the result of a least-mean-square fit to a Poisson distribu- 
tion, which is characterized by (Nt) = 3.0 (A = 2.8). The reduced average reflects the fact 
that the Poisson distribution underrepresents the tail of the distribution extracted from the 
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FIG. 6: Most frequent rupture force for single bonds as function of loading rate. The crosses 
are the maxima of the fitted single bond rupture force distributions D\{F) for the four different 
loading rates. The values are F* = 52.8, 64.4, 70.5 and 81.4 pN for increasing m. The dynamic 
force spectrum reveals a straight line with slope Fq = 7.85 pN which intersects F* = at fco-^b = 
0.159pN/s (marked by a square) so that ko = 0.02 Hz. The reactive compliance is Xj, = ksT/Fo = 
0.52 nm. 

experimental data. However, it is also important to note that the reconstruction at large 
forces is not very reliable due to the few events in the experimental data (compare fig. [3]). 
In order to test whether the Poisson distribution is compatible with our experimental data, 
we generated artificial rupture histograms by drawing random values for N t from the Pois- 
son distribution exactly as often as experimental data points were present. Fig. [5] shows a 
corresponding histogram for loading rate m = 1250 pN/s. The lines and the inset represent 
the same analysis as done before for the experimental data. The new Poisson distribution 
is slightly different, with A = 3.1 instead of A = 2.8, while the extracted values for fc an d 
F Q are essentially unchanged. Taken together, the data simulation suggests that the finite 
size of the data set can explain the deviations between experimental data and model fit. 

Figs.[3]and[4]suggest that the most salient feature of both the experimental and simulation 
data is the first peak, which corresponds to single bond rupture. In Fig. [6] we plot the fitted 
positions of the first peaks as a function of mm. As expected for dynamic force spectra, the 
data points can be fitted to a linear curve with the slope F = 7.85 pN and the intersection 
with F* = at k F = 0.159 pN/s. This again yields k = 0.02 Hz and Xb = 0.52 nm for 
the unstressed off-rate and the reactive compliance, respectively, demonstrating that our 
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analysis is consistent. These results suggest that the first peaks in a series of multiple bond 
rupture force histograms at different loading rates as shown in fig. [3] can be used to achieve 
a quick and simple estimate of the molecular parameters. In particular, this estimate can 
then be used as a starting point for the detailed numerical analysis presented above. 

V. DISCUSSION 

We have used multiple bond DFS-measurements with the BFP-setup to probe the outer 
barrier of the biotin-streptavidin bond. If multiple attachments are allowed, the frequency 
of successful binding events is larger than in classical single bond experiments so that in- 
formation can be obtained more efficiently in the form of rupture force histograms. As 
demonstrated here, analysis of these histograms requires theoretical modelling of the rup- 
ture of parallel adhesion bonds under an applied load because the peaks of the histograms 
follow from a superposition procedure. Using an established stochastic model for bond rup- 
ture under shared loading in combination with the assumption of a variable number of initial 
bonds, we showed for the first time how this superposition looks in practise, compare fig. HI 
The molecular bond parameters k Q and %h were extracted in good agreement with previous 
results from single bond experiments. 

One particular feature of our approach is that with a simple model for multiple bond 
rupture we arrived at an extraction of single bond data which is both successful and ef- 
ficient. Our model for force transduction and adhesion bonds is rather generic and can 
be easily extended to include additional effects like finite transducer and tether stiffnesses 
or generalised models for the binding landscape. E.g. it was shown for protein unfolding 
that analysis of histograms can be used to identify the shape of the binding landscape and 
identify deviations from the assumption of sharp transition barriers underlying DFS j^ . 
We expect that the agreement between experimental and theoretical histograms can be im- 
proved by including more substructure into the model. Recently, it has been proposed that 
deviations in rupture force histograms from the DFS predictions are caused by an intrinsic 
heterogeneity of bond parameters j^sj . The increased width and large force tails could be 
described imposing Gaussian noise on Fq. Although this effect will be small compared to 
the variations caused by the distribution of N t and will not alter the basic results, it might 
explain residual events at large rupture forces (another possible explanation is that some 
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initial population of the biotin-streptavidin ground state exists after all). Both examples 
show that for detailed studies of adhesion bonds, a model-based analysis of the histograms 
is required. 

In this work, rebinding effects could be safely neglected. It has been shown before that 
rebinding is only relevant if the dimensionless loading rate m/(koF ) < N t jiol JjJ]. Using 
the smallest value for the loading rate, m = 130pN/s, and the extracted bond parameters 
ko = 0.02 Hz and Fq = 7.83 pN, we get a value of 830, which is much larger than the largest 
value N™ = 10 used here. In principle, multiple bond experiments also offer the chance to 
measure rebinding rates, in marked contrast to single bond experiments. This could be done 
in different ways. For example, one could control the time allowed for bond formation prior 



to loading 



231 ] . The mean number of formed bonds could be extracted as described here and 



then converted into a value for the rebinding rate. Alternatively, one could work at small 
oading rate. Then the master equation from eq. ([1]) had to be extended by rebinding terms 
1^ | . A similar procedure as described here then could be used to extract the rebinding 
rate. For biotin-streptavidin, the above estimate suggests that then loading rates had to be 
smaller by at least one order of magnitude. Therefore adhesion bonds with faster off-kinetics 
than biotin-streptavidin might be more appropriate for this purpose. 

Finally, investigating multiple bond rupture will help to understand properties of biolog- 
ical adhesion sites which usually consist of clusters of adhesion molecules. For example, it 
has been found with image correlation microscopy that in living cells during cell migration, 
integrin adhesion receptors are preclustered with an average cluster size of three to four, 
which is very similar to the average cluster size studied here {2^. Thus biomimetic studies 
like the one presented here are essential to understand the way biological systems make use 
of adhesion clusters. 
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